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Abstract: AMANDA is a high-energy neutrino telescope presently under construction at the 
geographical South Pole. In the Antarctic summer 1995/96, an array of 80 optical modules 
(OMs) arranged on 4 strings (AMANDA-B4) was deployed at depths between 1.5 and 2 km. In 
this paper we describe the design and performance of the AMANDA-B4 prototype, based on 
data collected between February and November 1996. Monte Carlo simulations of the detector 
response to down-going atmospheric muon tracks show that the global behavior of the detector is 
understood. We describe the data analysis method and present first results on atmospheric muon 
reconstruction and separation of neutrino candidates. The AMANDA array was upgraded with 
216 OMs on 6 new strings in 1996/97 (AMANDA-BIO), and 122 additional OMs on 3 strings in 
1997/98. 
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1 Introduction 

Techniques are being developed by several groups to use high energy neutrinos as a probe for 
the highest energy phenomena observed in the Universe. Neutrinos yield information comple- 
mentary to that obtained from observations of high energy photons and charged particles since 
they interact only weakly and can reach the observer unobscured by intervening matter and 
undeflected by magnetic fields. 

The primary mission of large neutrino telescopes is to probe the Universe in a new observational 
window and to search for the sources of the highest energy phenomena. Presently suggested 
candidates for these sources are, for instance, Active Galactic Nuclei (AGN) and Gamma Ray 
Bursts (GRB). A neutrino signal from a certain object would constitute the clearest signature of 
the hadronic nature of that cosmic accelerator |Q]. Apart from that, neutrino telescopes search 
for neutrinos produced in annihilations of Weakly Interacting Massive Particles (WIMPs) which 
may have accumulated in the center of the Earth or in the Sun. WIMPS might contribute 
to the cold dark matter content of the Universe, their detection being of extreme importance 
for cosmology [Q]. Neutrino telescopes can be also used to monitor the Galaxy for supernova 
explosions y] and to search for exotic particles like magnetic monopoles Q. In coincidence with 
surface air shower arrays, deep neutrino detectors can be used to study the chemical composition 
of charged cosmic rays. Finally, environmental investigations - oceanology or limnology in water, 
glaciology in ice - have proved to be exciting applications of these devices JS], |^ . 

Planned high-energy neutrino telescopes differ in many aspects from existing underground neu- 
trino detectors. Their architecture is optimized to achieve a large detection area rather than 
a low energy threshold. They are deployed in transparent "open" media like water in oceans 
or lakes, or deep polar ice. This brings additional inherent technological challenges compared 
with the assembly of a detector in an accelerator tunnel or underground cavities. Neutrinos are 
inferred from the arrival times of Cherenkov light emitted by charged secondaries produced in 
neutrino interactions. The light is mapped by photomultiplier tubes (PMTs) spanning a coarse 
three-dimensional grid. 

The traditional approach to muon neutrino detection is the observation of upward moving muons 
produced in charged current interactions in the rock, water or ice below the detector. The 
Earth is used as a filter with respect to atmospheric muons. Still, suppression of downward- 
going muons is of top importance, since their flux exceeds that of upward-going muons from 
atmospheric neutrinos by several orders of magnitude. 

An array of PMTs can also be used to reconstruct the energy and location of isolated cascades 
due to neutrino interactions. Burst-like events, like the onset of a supernova, might be detected 
by measuring the increased count rates of all individual PMTs. 

Technologies for underwater telescopes have been pioneered by the since decommissioned DU- 
MAND project near Hawaii [^, 10| and by the Baikal collaboration |^, |8[. In contrast to these 



approaches, the AMANDA detector ||I^ uses deep polar ice as target and radiator. Two projects 



in the Mediterranean, NESTOR |11| and ANTARES [12], have joined the worldwide effort to- 
wards large-scale underwater telescopes. BAIKAL and AMANDA are presently taking data 
with first stage detectors. 



The present paper describes results obtained with the first four (out of the current thirteen) 
strings of the AMANDA detector. The paper is organized as follows: In section |2| we give a 
general overview of the AMANDA concept. Section ^ summarizes the results obtained with 
a shallow survey detector called AMANDA-A. Section ^ describes the design of the first four 
strings of the deeper array AMANDA-B4. Calibration of time response and of geometry are 
explained in section 5. In section we describe the simulation and reconstruction methods with 
respect to atmospheric muons and compare experimental data to Monte Carlo calculations. 
Section |^ demonstrates the performance of AMANDA-B4 operated in coincidence with SPASE, 
a surface air shower array. In section j^, the angular spectrum of atmospheric muons is derived 
and transformed into a dependence of the vertical intensity on depth. Section |9| describes the 
separation of first upward going muon candidates. Finally, a summary of the status of AMANDA 



and results is presented in section 10 



2 The AMANDA Concept 

AMANDA (Antarctic Muon And Neutrino Detector Array) uses the natural Antarctic ice as 
both target and Cherenkov medium. The detector consists of strings of optical modules (OMs) 
frozen in the 3 km thick ice sheet at the South Pole. An OM consists of a photomultiplier 
in a glass vessel. The strings are deployed into holes drilled with pressurized hot water. The 
water column in the hole then refreezes within 35-40 hours, fixing the string in its final position. 
In our basic design, each OM has its own cable supplying the high voltage (HV) as well as 
transmitting the anode signal. The components under the ice are kept as simple as possible, all 
the data acquisition electronics being housed in a building at the surface. The simplicity of the 
components under ice and the non-hierarchical structure make the detector highly reliable. 

Fig. 1^ shows the current configuration of the AMANDA detector. The shallow array, AMANDA- 
A, was deployed at a depth of 800 to 1000 m in 1993/94 in an exploratory phase of the project. 
Studies of the optical properties of the ice carried out with AMANDA-A showed that a high 
concentration of residual air bubbles remaining at these depths leads to strong scattering of light, 
making accurate track reconstruction impossible Q. Therefore, in the polar season 1995/96 a 
deeper array consisting of 80 OMs arranged on four strings (AMANDA-B4) was deployed at 
depths ranging from 1545 to 1978 meters, where the concentration of bubbles was predicted to 
be negligible according to extrapolation of AMANDA-A results. The detector was upgraded 
in 1996/97 with 216 additional OMs on 6 strings. This detector of 4+6 strings was named 
AMANDA-BIO and is sketched at the right side of fig. |. AMANDA-BIO was upgraded in the 
season 1997/98 by 3 strings instrumented between 1150 m and 2350 m which fulfill several tasks. 
Firstly, they explore the very deep and very shallow ice with respect to a future cube kilometer 
array. Secondly, they form one corner of AMANDA-II which is the next stage of AMANDA with 
altogether about 700 OMs. Thirdly, they have been used to test data transmission via optical 
fibers. 

There are several advantages that make the South Pole a unique site for a neutrino telescope: 

• The geographic location is unique: A detector located at the South Pole observes the 
northern hemisphere, and complements any other of the planned or existing detectors. 
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Figure 1: Scheme of the 1998 AMANDA instaUations. The left picture is drawn with true scahng. 
A zoomed view on AMANDA-A (top) and AMANDA-BIO (bottom) is shown at the center. The right 
zoom depicts the optical module. 



Ice is a sterile medium. The noise is given only by the PMT dark noise and by K"^^ decays 
in the glass housings, which are 0.5-1.5 kHz for the PMTs and spheres we used. Ocean 
and lake experiments have to cope with 100 kHz noise rates due to bioluminescence or 
K'^^ decays (25-30 kHz if normalized to the photocathode area of the 8" PMT used in 
AMANDA). This fact not only facilitates counting rate experiments like the search for 



low energy neutrinos from supernovae or GRBs, but also leads to fewer accidental hits in 
niuon events - an essential advantage for trigger formation and track reconstruction. 

• AMANDA can be operated in coincidence with air shower arrays located at the surface. 
Apart from complementing the information from the surface arrays by measurements of 
muons penetrating to AMANDA depths, the air shower information can be used to cali- 
brate AMANDA. 

• The South Pole station has an excellent infrastructure. Issues of vital importance to run 
big experiments like transportation, power supply, satellite communication and technical 
support are solved and tested during many years of operation. Part of an existing building 
can be used to house the surface electronics. 

• The drilling and deployment procedures are tested and well under control. AMANDA 
benefits from the drilling expertise of the Polar Ice Coring Office (PICO). Currently about 
five days are needed to drill a hole and to deploy a string with PMTs to a depth of 2000 m. 
Future upgrades of the drilling equipment are expected to result in a further speed-up. 

The optical properties of the ice turned out to be very different from what had been expected 
before the AMANDA-A phase. Whereas absorption is much weaker than in oceans, scatter- 
ing effects turned out to be much stronger. Even at depths below 1400 meters, where residual 
bubbles have collapsed almost completely into air hydrates, scattering is nearly an order of mag- 
nitude stronger than in water (see below) . Since scattering of light smears out the arrival times 
of Cherenkov flashes, a main question was whether under these conditions track reconstruction 
was possible. As shown below, the answer is yes. 



3 AMANDA-A: A First Survey 

Preliminary explorations of the site and the drilling technology were performed in the Antarc- 



tic Summer 1991/92 |13]. During the 1993/94 campaign, four strings each carrying 20 OMs 
("AMANDA-A") were deployed between 800 and 1000 m depth. None of the 73 OMs (equipped 
with 8" EMI PMTs) surviving the refreezing process failed during the following two years, giv- 
ing a mean time between failures (MTBF) > 40 years for individual OMs in AMANDA-A. The 
OMs are connected to the surface electronics by coaxial cables. Along with the coaxial cables, 
optical fibers carry light from a Nd:YAG laser at the surface to nylon light diffusers placed 
about 30 cm below each PMT (see fig. ||). Time calibration is performed by sending nanosecond 
laser pulses to individual diffusers and measuring the photon arrival time distributions at the 
closest PMT. From the distribution of the arrival times at distant PMTs, the optical properties 
of the medium were derived |^, ^ . The measured timing distributions indicated that photons do 
not propagate along straight paths but are scattered and considerably delayed due to residual 
bubbles in the ice. The distributions could be fitted well with an analytic function describing 
the three-dimensional random walk (scattering) including absorption. These results showed that 
polar ice at these depths has a very large absorption length, exceeding 200 m at a wavelength of 
410 nm. Scattering is described by the effective scattering length L^fj = Lgc/iX ~ (cos 0)), where 
Lsc is the geometrical scattering length and {cos 6) the average cosine of the scattering angle Q. 



Le// increases with depth, from 40 cm at 830 m depth to 80 cm at 970 m. In accordance with 



measurements at the Vostok Station (East Antarctica |14|) and Byrd Station (West Antarctica) 
these results suggested that at depths greater than 1300-1400 m the phase transformation from 
bubbles into air-hydrate crystals would be complete and bubbles would disappear. 

Although not suitable for track reconstruction, AMANDA-A can be used as a calorimeter for 



energy measurements of neutrino-induced cascade-like events |1(:]. It is also used as a supernova 



monitor [15|. Events that simultaneously trigger AMANDA-A and the deeper AMANDA-B 
have been used for methodical studies like the investigation of the optical properties of the ice 
or the assessment of events with a lever arm of one kilometer. 



4 Deployment and Design of AMANDA-B4 

4.1 Drilling and Deployment Procedure 

Drilling is performed by melting the ice with pressurized water at 75°C. The drilling equipment 
operates at a power of 1.9 MW and the typical drill speed is about 1 cm/s. It takes about 3.5 
days to drill a 50-60 cm diameter hole to 2000 m depth. 

In the season 1995/96, we drilled four holes, the deepest of them reaching 2180 m. It took 
typically 8 hours to remove the drill and the water recycling pump from the completed hole. 
The deployment of one string with 20 OMs and several calibration devices took about 18 hours 
(with a limit of 35 hours set by the refreezing of the water in the hole) . 

Several diagnostic devices allow monitoring of the mechanical and thermal parameters during 
the entire refreezing process and afterwards. It was shown that the temperature increases with 
depth in good agreement with the prediction of a standard heat flow calculation for South Pole 
ice. At the greatest depth, the temperature of the ice is ~ -31°C, about 20° warmer than at the 
surface. During the refreezing, the pressure reached a maximum of 460 atm, more than twice 
the hydrostatic pressure which is asymptotically established. 

4.2 Detector Design 

The four strings of AMANDA-B4 were deployed at depths between 1545 and 1978 m. An OM 
consists of a 30 cm diameter glass sphere equipped with a 8" Hamamatsu R5912-2 photomulti- 
plier, a 14-dynode version of the standard 12-dynode R5912 tube. The PMTs are operated at 
a gain of 10^ in order to drive the pulses through 2 km of coaxial cable without in-situ amplifi- 
cation. The amplitude of a one-photoelectron pulse is about 1 V. The coaxial cable is also used 
for the HV supply, with the advantage that only one cable and one electrical penetrator into 
the sphere are required for each OM. The measured noise rate of the AMANDA-B4 PMTs is 
typically 400 Hz (threshold 0.4 photoelectrons). 

The photocathode is in optical contact with the glass sphere by the use of silicon gel. The 
transmission of the glass of the pressure sphere is about 90% in the spectral range between 400 
and 600 nm; the 50% cutoff on the UV side is at about 365 nm. The glass spheres are designed 
to withstand pressures of about 660 atm. 
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Figure 2: AMANDA-B4: Top view, witli distances between strings given in meters, and side view 
showing optical modules and calibration light sources. Upward looking PMTs are marked by arrows. 



Each string carries 20 OMs with a vertical spacing of 20 m. The fourth string carries six addi- 
tional OMs connected by a twisted pair cable. These six OMs will not be used in the analyses 
presented in this paper. 



Fig. ^ shows a schematic view of AMANDA-B4. Ah PMTs look down with the exception of jj^ 
1,10 in strings 1 to 3 and #1,2,10,19,20 in string 4 (with the numbers running from top to bottom 
of a string). Strings 1-3 form a triangle with side lengths 77-67-61 m; string 4 is close to the 
center. The OMs are arranged at depths 1545-1925 m (string 1), 1546-1926 m (string 2), 1598- 
1978 m (string 3) and 1576-1956 m (string 4). The additional six OMs equipped with twisted pair 
cables are at string 4 between 2009 and 2035 m. Seven of the 80 PMTs which define AMANDA- 
B4 were lost due to overpressure and shearing forces to the electrical connectors during the 
refreezing period. These losses can be reduced by computer controlled drilling avoiding strong 
irregularities in the hole diameter, and by improved connectors. Another 3 PMTs failed in the 
course of the first 3 years of operation, giving a MTBF of 73 years. 

4.3 Electronics and DAQ 

Each PMT can give a series of pulses which can be resolved if separated from each other by more 
than a few hundred nanoseconds. The data recorded consist of the leading and trailing edges 
of the pulses. The time-over-threshold gives a measure of the amplitude of individual pulses. 
Another measure of the amplitude is obtained by a voltage sensitive ADC which records the 
peak value out of the subsequent hits of an event in a PMT. Actually, the information consists 
of leading and trailing edges of the last 8 resolved pulses, and of the largest amplitude of those 
of them which lie in a 4 ^usec window centered at the array trigger time. Also recorded is the 
GPS time at which the event occurred. A scheme of the AMANDA electronics layout is shown 
in fig. H 
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Figure 3: DAQ system used for AMANDA-B4 during 1996 



The signal from each cable is fed to a module consisting of a DC blocking high-pass filter which 
picks up the pulse, a fan-out sending it to 2 amplifiers with 100 x and 25 x gain, and a 2 /xsec 
delay for the low-gain signal. 

The delayed signal is sent to a Phillips 7164 peak sensing ADC. The other pulse is split and sent 
to LeCroy 4413 discriminators with thresholds set at 100 mV corresponding to about 0.3-0.4 
photoelectrons at the given high voltage. One of the resulting ECL pulses is fed into a LeCroy 
3377 TDC while the other is sent to the majority trigger. The TDC records the last 16 time 
edges occurring within a 32 ^usec time window. 

The majority logic requests > 8 hit PMTs within a sliding window of 2 /xsec. The trigger 
produced by this majority scheme is sent to the NIM trigger logic. The latter accepts also 
triggers from AMANDA-A or the air shower experiments SPASE-1, SPASE-2 and GASP. Thus 
AMANDA also records data when these detectors trigger even if a proper AMANDA trigger is 
not fulfilled. The total trigger rate during 1996 was about 26 Hz on average. The coincidences 
from the other detectors contributed about 8 Hz to the total rate. 

The differences in cable length are not compensated before triggering. Therefore the true trigger 
window would be about 300 nsec for a vertically downgoing relativistic particle and ~ 4/isec for 
an upgoing one. As a result downgoing particles are suppressed compared to upgoing. 

Upon triggering, an ADC gate of 4 /isec width is formed, a stop signal is sent to the TDCs and a 
readout signal is sent to a Hytec LP1341 list processor. Then a veto lasting several microseconds 
inhibits further trigger signals. 

A separate system ("SN scalers" in fig. y) monitors the counting rates of individual PMTs 
and searches for rate excesses lasting several seconds. Such an increase would be expected for 
multiple low-energy neutrino interactions close to each PMT due to a supernova burst |^, |l^ ]. 

The AMANDA-B4 DAQ was running on a Macintosh Power PC communicating through a SCSI 
bus with the CAMAC crate controller. From the distribution of the time differences between 
subsequent events, the dead time of the DAQ is estimated to be about 12%. The Macintosh 
has been replaced by a Pentium-II PC running under LINUX in 1998, and part of the CAMAC 
electronics by VME modules. 

Fig. shows the distribution of the leading-edge times of one PMT for data taken with the 
8-fold majority trigger. The sharp peak at 23 //sec is given by the time when this PMT was the 
triggering one (i.e. the eighth) within a 2/isec window. The fiat part is due to noise hits and 
the bulge after the main distribution to afterpulses (about 6%.) 



4.4 Calibration Light Sources and Ice Properties 

An essential ingredient to the operation of a detector like AMANDA is the knowledge of the 
optical properties of the ice, as well as a precise time calibration of the detector. Various light 
calibration sources have been deployed at different depths in order to tackle these questions: 
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Figure 4: Leading edge times of PMT # 10 of AMANDA-B4 for data taken witlr an 8-fold majority 
trigger. 

• The YAG laser calibration system. It uses optical fibers with difi'users located at each 
PMT. This system is similar to that used for AMANDA-A. The range of transmittable 
wavelengths is > 450 nm, the time resolution is about 15nsec at 530 nm, the maximum 
intensity emitted by the diffusers is 10^ photons/pulse. Apart from ice investigations, the 
laser system is used for time calibration of the PMT closest to the diffuser and for position 
calibration (see section ||). 

• A nitrogen laser at 1850 m depth, wavelength 337 nm, pulse duration 1 nsec, with a 
maximum intensity of 10^^ photons/pulse. 

• Three DC halogen lamps (one broadband and two with filters for 350 and 380 nm), 
maximum intensity lO^'' (UV-filtered) and 10^^ (broadband) photons/second. 

• LED beacons, operated in pulsed mode (500 Hz, pulse duration 7 nsec, 10^ photons/pulse) 
and DC mode (10^"^ to 10^^ photons/sec), wavelength 450 nm. A filter restricts the output 
of a few beacons to 390 nm, with reduced intensity. 



Time-of-flight measurements have been made for a large variety of combinations of optical fiber 
emitters and PMTs for the YAG laser system, and at different wavelengths and intensities. The 
nitrogen laser provided data at 337 nm. The result is a considerable data base of hundreds of 
time distributions. The width of the distributions is sensitive predominantly to scattering and 



the tail to absorption (see [29| for details). The DC sources provide data for attenuation, i.e. 
the combined effect of absorption and scattering. 
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The YAG laser results indicate a dramatic improvement compared to AMANDA- A results. 
Fig. |5| shows the distributions of arrival time for source-detector distances of 20 and 40 m, 
respectively, for AMANDA-A as well as AMANDA-B depths. The much smaller widths for 
AMANDA-B support the expectation that bubbles as the dominant source of scattering have 
mostly disappeared at depths between 1550 and 1900 m IIT4 
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Figure 5: Arrival time distributions for 510 nm photons for two source-detector distances. Black his- 
tograms: AMANDA-B. Hatched histograms:AMANDA-A. The histograms are normalized to the same 
area. 



Details of the analysis of the optical properties of the ice at AMANDA-B4 depths have been 
published elsewhere [0]. Final results will be published in a separate paper. Figure y shows 
preliminary data on the wavelength dependence of the coefficients for scattering, be, and ab- 
sorption, a. The absorption length \a = 1/a is between 90 and 100 m for wavelengths below 
460 nm, i.e. ice does not degrade in transparency towards smaller wavelengths down to 337 
nm. The effective scattering length Ag// = l/&e varies between 24 and 30 m in the relevant 
wavelength range. Ag// = Xscatt/{^ — {cos9)), with Xscatt being the geometric scattering length. 
(cos 6) is the average cosine of the scattering angle and is supposed to be about 0.8 in deep ice. 
The attenutation length Xatt which characterizes the decrease of the photon flux as a function 
of the distance is about 27 m. These values are averages over the full depth interval covered by 
AMANDA-B4. The variation of attenuation over this depth range is within ±30%. 



11 




600 



wavelength (nm) 



Figure 6: Absorption (a) and scattering (bg) coefficients at an average depth of 1.7 km, compared to 
theory of He and Price |jl3 . 

5 Calibration of Time Response and Geometry 

5.1 Time Calibration 

The measured arrival times from each PMT have to be corrected for the time offset to, that is, 
the time it takes a signal to propagate through the PMT and the coaxial cable and get digitized 
by the DAQ. The time offset is determined by sending light pulses from the YAG laser to the 
diffuser nylon balls located below each OM. Two fibers are available for each PMT, one single 
and one multi-modal. The time it takes for light to travel though the fiber is measured using 
an OTDR (Optical Time Domain Reflectometer) and subtracted from the time distributions 
recorded. 

For each PMT, the time difference between the laser pulse at the surface and the PMT response 
arriving back is measured. Upon arrival at the surface, the pulses have traveled through nearly 
2000 meters of cable and are dispersed, with typical time-over-thresholds of 550 nsec and rise 
times of 180 nsec. The threshold used for TDC measurements is set to a constant value with 
the consequence that small pulses will reach that value later than larger ones. This causes an 
amplitude-dependent offset or "time walk", which can be corrected for by 



ttrue = tlE -to- a/VADC. 



(1) 
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Here, Ile is the measured leading edge time and ttme the true time at which the hght pulse 
reaches the photocathode. The estimates of the time offset io and the time-walk term a are 
extracted from scatterplots like the one shown in fig. |^ 
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Figure 7: Example of a fitted leading edge (with 100<ADC<1200) for module 19 on string 3. The ADC 
value measures the peak value of the amplitude. 



The time resolution achieved in this way can be estimated by the standard deviation of a 
Gaussian fit to the distribution of time residuals after correction, yielding 4-7 nsec (see Fig. 1^ 
for an OM with 4 nsec resolution). Part of the variation is due to quality variations of the 1996 
optical fibers. Laboratory measurements yield a Gaussian width of 3.5 nsec after 2 km cable. 
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Figure 8: Residuals after subtracting the time correction obtained with the fitted parameters io and a 
for module 19 on string 3. The standard deviation of the Gaussian fit is 4 nsec. 
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5.2 Position Calibration 

Information about the exact geometry of the array can be obtained by different methods. Firstly, 
the measured propagation times of photons between different hght emitters and receivers can 
be used to determine their relative positions. Secondly, absolute positions can be obtained from 
drill recordings and pressure sensors. 

Laser Calibration 

The YAG laser, the nitrogen laser and the pulsed LEDs can be used to infer the OM positions 
from the time-of-flight of photons between these light sources and the OMs. The zero time 
is determined from the response of the OM closest to the light source which is triggered by 
unscattered photons. This PMT is lowered in voltage in order not to be driven in saturation, 
and a time correction accounting for the longer PMT transit time is added. In contrast to the 
close OM, the distant OMs see mostly scattered photons. However, for a few of the events out 
of a series of about 1500 laser pulses, the leading edge should be produced by photons which are 
only slightly scattered. Therefore the distance between emitter and OM can be estimated from 
the earliest events in the time-difference distribution (see fig. ^). 
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Figure 9: Simulated time-shift distribution for 1500 one-photoelectron events, for a distance of 60 m 
between emitter and receiver. A Gaussian smearing of lOnsec was applied to individual entries. Clear 
ice would yield a lOnsec wide peak at Onsec. 



In order to reduce the sensitivity to fluctuations in the number of early hits and binning effects, 
the whole left flank of the distribution is fltted with a Gaussian between the maximum of the 
distribution (heightO in flg. H) and the flrst bin with a height larger than height 1 =1/10 
heightO. The corrected "flrst" time is given by that bin (binl) for which the fltted Gaussian 
yields a height exceeding heightl. This time has to be corrected further for the shifts due 
to scattering which are expected even for the flrst bin of the distribution. The corrections 
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were obtained from Monte Carlo (MC) calculations and are almost insensitive to variations in 
absorption and scattering length of a few meters. 

Given the limited number of measured emitter-OM combinations available for AMANDA-B4, it 
would have been impossible to keep the coordinates of each OM as free parameters in a global 
position fit. Therefore, all strings were assumed to be straight and parallel and the OMs to be at 
a fixed vertical distance (20 m) relative to each other. For each emitter covering enough OMs, a 



graph of the distance d{zi) between source and OM i versus depth Zi can be drawn (see fig. 10). 
The inter-string distance D and emitter depth zo with respect to the Zi can be estimated from 
this graph by fitting (fig. [ill) : 



d{zi) = ^D^ + {z, - zo)2. (2) 

The residuals from all fits to the 1996 AMANDA-B4 data have a standard deviation of 2 m. 




^D'HzrZo)' 



receivers 



Figure 10: Principle of position measurement 



In 1996-1997, six more strings were added on the outside of the B4 detector, and a new position 
calibration performed. The increased statistics and possibilities of new cross-checks and con- 
straints enabled correction of the existing geometry with an uncertainty of 1 m in the horizontal 
plane and 0.5-1.0 m in depth. 

Drill data 

The geometry of the array is surveyed in an independent way by monitoring the position of the 
drill-head while it is going down each hole. The data were recorded by the drill instrumentation 
at each 10 cm step, recording the path-length, the value of the Earth's magnetic field as measured 
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Figure 11: Fit of the distance d{zi) versus depth-shift Zi — zi between OMs at string 4 and a laser 
emitter at string 2. String distance D and depth shift zq — Z2 are given by the minimum of the parabola. 



by a flux magnetometer and the angles (bank and elevation) given by perpendicular pendulums. 
This information can then be used to reconstruct the hole profiles. The results found are 
compatible with the laser measurements within 1-2 m in the horizontal plane. The advantage 
of this method is that it yields positions relative to the surface, i.e. in a global reference frame. 
It also takes into account tilts in the strings. However, it does not yield the depth locations of 
the OMs. The absolute depths of the strings were given by pressure sensors deployed with the 
OMs. 



6 Simulation and Reconstruction of Muons 



6.1 Simulation 



Downgoing muons are generated by full atmospheric shower programs which simulate the pro- 
duction of muons by isotropic primary protons [19| or protons and nuclei [^] with energies up 
to 1 PeV. The muons are propagated down to a plane close to the detector. Upgoing muons 



are generated from atmospheric neutrinos, using the flux parameterization given in |21], from 
neutralinos annihilating in the center of the Earth, using the flux calculations of [^, and from 
point sources, using arbitrary energy distributions and source angles; they may start anywhere 
within the fiducial volume (which increases with increasing neutrino energy due to the muon 
range) and are propagated simulating the full stochastic energy loss according to |22]. 
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It would be computationally impractical to generate and follow the path of each of the multiply 
scattered Cherenkov photons produced by muons and secondary cascades for every simulated 
event. Therefore, this step is accomplished by doing the photon propagation only once by 
a separate MC program and storing the results in large multidimensional tables. The tables 
give the distribution of the mean number of photoelectrons expected and of the time delay 
distribution, as a function of the position and the orientation of a PMT relative to the muon 
track. They include the effects of the wavelength dependent quantum efficiency, the transmission 
coefficients of glass spheres and optical gel, and the absorption and scattering properties of the 
ice. Once the tables are compiled, events can be simulated quickly by locating the PMT relative 
to any input particle and looking up the expected number and time distribution of photoelectrons 
in the tablesjj. The known characteristics of the AMANDA PMTs, the measured pulse shapes, 
pulse heights and delays after signal propagation along the cables, and the effect of electronics 



are then used to generate amplitude and time information [23| 



6.2 Reconstruction 

The reconstruction procedure for a muon track consists of five steps: 

1. Rejection of noise hits, i.e. hits which have either a very small ADC value or which are 
isolated in time with respect to the trigger time or with respect to the nearest hit OM. 



2. A line approximation following |24] which yields a point on the track, r, and a velocity v: 

{nti) - {ri){ti) 



r= (ri) -V- (ti) 



(tl) - {Uf 



with fi and ij being the coordinate vector and response time of the i-th PMT. 

3. A likelihood fit based on the measured times which takes the the track parameters obtained 
from the line fit as start values. This "time fit" yields angles and coordinates of the track 
as well as a likelihood Cume- 

4. A likelihood fit using the fitted track parameters from the time fit and varying the light 
emission per unit length until the probabilities of the hit PMTs to be hit and non-hit 
PMTs to be not hit are maximized. This fit does not vary the direction of the track but 
yields a likelihood Chu with can be used as a quality parameter. 

5. A quality analysis, i.e. application of cuts in order to reject badly reconstructed events. 
Steps 3 and 5 are outlined in the following two subsections. 

6.3 Time Fit 

In an ideal medium without scattering, one would reconstruct the path of minimum ionizing 
muons most efficiently by a x^ minimization process. Due to scattering in ice, the distribution 

^This method assumes that ice is isotropic and and homogeneous which is reasonable in a first approximation: 
firstly, since the variations of the original ice with depth have been measured to be smaller than ±25%, secondly, 
since the freshly frozen ice in the holes occupies only a small volume of the array. 
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Figure 12: Delay-time distributions for modules facing (full curves) and back-facing (dashed curves) a 
muon track. Data are shown for muon tracks with impact parameters of 5 meters (a) and 150 meters 
(b). 



of arrival times of photoelectrons seen by a PMT is not Gaussian but has a long tail at the high 



side - see fig. 12. 



To cope with the non-Gaussian timing distributions we used a likelihood analysis. In this 
approach, a normalized probability distribution function pi(t) gives the probability of a certain 
time delay t for a given hit i with respect to straightly propagating photons. This probability 
function is derived from the MC simulations based on the photon propagation tables introduced 



in section 3.1. The probability depends on the distance and the orientation of the PMT with 
respect to the muon track. By varying the track parameters the logarithm of a likelihood 
function C is maximized. 



18 



log{C) = log Yl pi\ = J2_ iog(^*) 

Vail hits / all hits 

In order to be used in the iteration process, the time delays as obtained from the separate 
photon propagation Monte-Carlo have to be parameterized by an analytic formula. The pa- 
rameterization of the propagation model itself is extended to allow for timing errors of PMTs 
and electronics as well as the probability of noise hits at random times. The AMANDA collab- 
oration has developed two independent reconstruction programs, which are based on different 
parameterizations of the photon propagation and different minimization methods p^, 28, 26 1. 



The comparison of these algorithms and the use of different optical models show that the results 
of both methods are in good agreement with each other and do not depend on a fine-tuning of 



the assumed optical parameters. Fig. |12| shows the result of the parameterization of the time 
delay for two distances and for two angles between the PMT axis and the muon direction. At a 
distance of 5 m and a PMT facing toward the muon track, the delay curve is dominated by the 
time jitter of the PMT. However, if the PMT looks in the opposite direction, the contribution 
of scattered photons yields a long tail towards large delays. At distances as large as 150 m, 
distributions for both directions of the PMT are close to each other since all photons reaching 
the PMT are multiply scattered. 

The parameterization used for most of the results presented in this paper is a Gamma distribution 
modified with an absorption term |27| 



^-{d/X) Ad/X-1) 

with the distance r between OM and muon track, the scaled distance d ~ 0.8/ sm{6c) ■ r, the 
absorption length Xq and only two parameters r ~ 700 ns and A ~ 50 m. 

The second approach uses an F-distribution with an exponential tail for large time-delays, which 
results in a comparable accuracy |2| 



6.4 Quality Analysis 

Quality criteria are applied in order to select events which are "well" reconstructed. The criteria 
define cuts on topological event parameters and observables derived from the reconstruction. 
Below we list those used in the following: 

• Speed \v\ of the line fit. Values close to the speed of light indicate a reasonable pattern of 
the measured times, values smaller than 0.1 m/nsec indicate an obscure time pattern. 

• "Time" likehhood per hit PMT log{£tvme)/Nhit- 

• Summed hit probability for all hit PMTs J2 Phit- 

• "Hit" likelihood normalized to all working channels, \og{Chit)/Nau- 

19 



The latter two parameters are good indicators of whether the location of the fitted track, 
which relies exclusively on the time information, is compatible with the location of the 
hits and non-hits within the detector. 

Number of direct hits, N^ir, which is defined to be the number of hits with time residuals 
ij (measured) — tj(fit) smaller than a certain cut value. We use cut values of 15 nsec, 25 nsec 
and 75 nsec, and denote the corresponding parameters as Ndir{^^), -^dir(25) and Ndiri75), 
respectively. Increasing the time window leads to higher cut values in N^ir but allows a 
finer gradation of the cut. 

Events with more than a certain minimum number of direct hits (i.e. only slightly delayed 
photons) are likely to be well reconstructed. This cut turned out to be the most powerful 
cut of all ||. 

The projected length of direct hits onto the reconstructed track, L^ir- A cut in this 
parameter rejects events with a small lever arm. 

Vertical coordinate of the center of gravity, zcoG- Cuts on this parameter are used to 
reject events close to the borders of the array. Very distant tracks are not likely to be well 
reconstructed. 



Fig. 13 shows the distribution of two of these observables, the number of direct hits within 
15 nsec, Ndir{^5), and the summed hit probability ^Phu of all hit channels. It demonstrates 
the good agreement between results from MC and experiment. 
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Figure 13: Distributions of two reconstructed event observables for MC down-going muon events (dashed 
lines) and from experimental data (full lines). Left: Number of direct hits, NfHr{\h)] Right: smnmcd hit 
probability, Y.Phit, 
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Fig. 14 demonstrates the effect of cuts on the number of direct hits and the summed hit proba- 
bility on the reconstructed angular distribution of experimental data and the MC sample. The 
cuts are Ndirecti^^) ^ 5 and J2 Phit ^ 2.5. Both samples are dominantly due to down-going 
atmospheric muons. Despite that, a small but similar fraction of events is falsely reconstructed 
as up-going events. After application of the above quality criteria the tail below the horizon 
almost disappears. Note that not only the shapes but also the absolute passing rates on all cut 
levels are in good agreement between data and Monte Carlo. The angular mismatch between the 
reconstructed muon angle and the original angle used in the MC simulation after both cuts is 5.5 
degrees. We note that this value strongly depends on the particular set of cuts, the minimum 
acceptable passing rate, the incident angle of the muon, and the range of muons stopping in the 
array. 
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Figure 14: Reconstructed zenith angle distributions of experimental data (line) and downward muon 
MC events (points) after a stepwise application of quality cuts. 
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7 SPASE coincidences 



AMANDA is unique in that it can be calibrated by muons with known zenith and azimuth 
angles which are tagged by air shower detectors at the surface. AMANDA-B4 has been running 
in coincidence with the two SPASE (South Pole Air Shower Experiment) arrays, SPASE-1 [^^ 
and SPASE-2 @. SPASE-1 was located 840 m from the center of the AMANDA array projected 
to the surface, whereas SPASE-2 is located 370 m away (see fig jig). The scintillation detectors 
of SPASE-2 are complemented by an array of air Cherenkov detectors pO| , ^. The primary 
goal of these devices is the investigation of the chemical composition of primary cosmic rays in 
the region of the "knee" [|32|| . Another detector, the gamma imaging telescope GASP, is also 
operated in coincidence with AMANDA. 
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Figure 15: Side view of the two SPASE arrays relative to AMANDA- A and AMANDA-B. 



In this section, we summarize calibration results obtained from the coincident operation of 
AMANDA and SPASE-2. SPASE-2 consists of 30 scintiUator stations of 0.8 m^ on a 30 m 
triangular grid. The area of the array is 1.6 • 10^ m'^, and it has been running since January 
1996. For each air shower, the direction, core location, shower size and GPS time are determined. 
Showers with sufficient energy to trigger SPASE-2 (~ 100 TeV) yield on average 1.2 muons 
penetrating to the depth of AMANDA-B. On every SPASE-1 or SPASE-2 trigger, a signal is 
sent to trigger AMANDA. The GPS times of the separate events are compared offiine to match 
coincident events. 
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A one- week sample of these events has been analyzed m order to compare the directions of muons 
determined by AMANDA-B4 to those of the showers measured by SPASE-2. A histogram of 



the zenith mismatch angle between SPASE-2 and AMANDA-B4 is shown in fig. 16. The selected 



events are required to have >8 hits along 3 strings and to yield a track which is closer than 
150m to the air shower axis measured by SPASE-2 (upper histogram). The hatched histogram 
shows the distribution of the zenith mismatch angle after application of the following quality 
cuts: 

• likelihood \og{Ctime) / Nhu > -12, 

• more than four hits with residuals smaller than 75 nsec {N^ir{75) > 4), 

• length of the projection of OMs with direct hits to the track larger than 50 meters 
{Ldir{75)>50m). 
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Figure 16: Mismatch between zenith angles determined in AMANDA-B4 and SPASE-2. 

428 of the originally 840 selected events pass these quality cuts. The Gaussian fit has a mean 
of (0.14 lb 0.19) degrees and a width of cr = (3.6 ± 0.17) degrees. This is nearly 2 degrees better 
than the resolution obtained in the previous section for all downward muons and for a different 
set of cuts. MC yields a resolution of about 4 degrees. 

The small mean implies that there is little systematic error in zenith angle reconstruction. The 
SPASE-2 pointing accuracy, which contributes to the average mismatch, depends on zenith angle 
and shower size. For most of the coincidence events, the SPASE-2 pointing resolution, defined 
as the angular distance within which 63% of events are contained, is between 1° and 2° [g^, 31|. 
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8 Intensity-vs-Depth Relation for Atmospheric Muons 

8.1 Angular Dependence of the Muon Flux 

In section g, the muon angular distribution was shown as a function of various cuts in order to 
demonstrate the agreement between experimental data and MC simulations. In this section, we 
calculate the muon intensity I as a function of the zenith angle 9. I{9^) is given by 

where 

• N^{9) is the number of muons assigned by the analysis to a zenith angle interval cen- 
tered around cos9^. For the analysis presented in this section, we start from the angular 
distribution Ni^i{6rec) obtained from the reconstruction, without applying cuts. This dis- 
tribution is strongly smeared (see fig.|lj, top). We have calculated the elements of the 
parent angular distribution Nfj_{6) from the reconstructed distribution N^{6rec) using a 
regularized deconvolution procedure [^, |3^ . 

• r is the run time. We used the data from June 24, 1996, with T = 22.03 hours, and 
9.86 • 10^ events triggering AMANDA-B4. 

• Sdead corrects for the dead time of the data acquisition system. This factor was determined 
from the time difference distribution of subsequent events. The dead-time losses for the 
two runs used in this analysis are 12%, i.e. Sdead = 1/0.88 = 1.14. 

• AQ is the solid angle covered by the corresponding cos 9^ interval. 

• ^e//(cut, 9fj_) is the effective area, after the application of a multiplicity trigger, for a given 
cut at zenith angle 0^. The effective area is shown in fig. ^ as a function of the zenith 
angle and for different cuts on the number of hit OMs. 

• erecidfi) is the reconstruction efficiency for zenith angle 9^ which ranges between 0.82 at 
cos 9 = 1.0 and 0.75 at cos 9 = 0.2. 

• m{9y) is the mean muon multiplicity at angle 9^ at the "trigger depth". The trigger depth 
heff was defined as depth of zqm-, the center of gravity in the vertical coordinate z of all 
hit OMs. The average /ig// depends on the angle. It is highest for cos9 between 0.4 and 
0.8 (about 30 m below the detector center) and falls toward the vertical (at maximum 
80m below the center). The mean muon multiplicity is about 1.2 for vertical tracks and 
decreases towards the horizon. Since the generator used in this analysis [^] simulates only 
proton induced showers, this value is an underestimation by about 10%. 

Fig. ^ shows the angular distribution of the flux of downgoing muons, I{9^), as obtained from 
eq.3. In order to illustrate the stability of the method with respect to cuts biasing the measured 
angular distribution, the flux is shown for samples defined by different majority triggers {Nhn > 
8, 10, 12, 14, 16, 18). Apart from the point closest to the horizon which is not only most strongly 
biased but also has the lowest statistics, deviations are within 25%. For further studies we use 
the sample with Nhu > 16. 
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Figure 17: Effective trigger area of AMANDA-B4 as a function of zenith angle, for 3 different majority 
criteria. 
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Figure 18: Angular distribution of the downward going muon flux, I{9^), as obtained from eq.(3). 
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8.2 Transformation of Angular Flux to Vertical Intensity as a Function of 
Depth 

The measured flux I {9) can be transformed into a vertical flux 1(6 = 0,h), where h is the ice 
thickness seen under an angle 9: 



l{e = 0,h) =I{9)-cos{9) 



(4) 



The cos ^-conversion correcting for the sec{9) behavior of the muon flux is valid for angles up to 
60° [46|. The term Ccorr taken from |Q| corrects for larger angles and lies between 0.8 and 1.0 
for the angular and energy ranges considered here. 



The vertical intensities obtained in this way are plotted in fig. 19 and compared to the depth- 
intensity data published by DUMAND |45] and Baikal ]|], and to the prediction by Bugaev et 



al. [38|. One observes satisfying agreement of all experiments with the prediction. 
We also fitted our data to a parameterization taken from 



,40]: 



I{h) =h-E^ 



crit 



'0 • 



^eff 



A^effh) 



(5) 



Ecrit is the minimum energy necessary to reach the depth h. It is obtained from the pa- 
rameterization flq] dE/dx = a + b ■ E^ where a ~ 2 MeV/(g-cm~^) denotes the continuous 
energy loss due to ionization, and b{E^) is proportional to the stochastic energy loss due to 
pair production, bremsstrahlung and nuclear cascades. Prom this parameterization one obtains 
Ecrit = 0,/b ■ [exp(5 ■ h) — 1]. Iq is the normalization parameter and 7 ~ 2.78 |4C] the spectral 
index. We approximate b{E^) by an energy independent parameter heff- Fitted to equation (5), 
our data for the vertical intensity result in the following values for Iq and b^ff'- 



Iq = (5.04 ±0.13) cm-^g- 



ster 



5e// 



(2.94 ±0.09) • 10-^g 



^^ cm^. 



This compares to Iq = (5.01±0.01) cm-^g-ister"^ and h^jf = (3.08±0.06) 10"^ g-^ cm^ obtained 
for Nhit > 8, showing that the result is rather insensitive to the actual cut condition. 

For the purpose of completeness we give also the results for the more usual parameterization 



Iih,9^ = 0) 



-h/X 



(6) 



where a is set to [41|, to 2 [42| or is a free parameter |43|. The purely exponential dependence 
(a = 0) clearly does not describe the data at depths smaller than 4-5 km. Leaving all parameters 
free e|, one obtains a,, = (0.89 ± 0.30) • 10"^ cm-2s-ister-\ A = (1453 ± 612) gcm'^^ and 



a = 2.0 ± 0.25, being also in agreement with a fixed as in [42|. 
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Figure 19: Vertical intensity versus depth for AMANDA, BAIKAL and DUMAND. The soUd hne gives 
the prediction of p8[ which coincides with the curves obtained from the parameterizations (5) and (6). 

9 Search for Upward Going Muons 

AMANDA-B4 was not intended to be a full-fledged neutrino detector, but instead a device 
which demonstrates the feasibility of muon track reconstruction in Antarctic ice. The limited 
number of optical modules and the small lever arms in all but the vertical direction complicate 
the rejection of fake events. In this section we demonstrate that in spite of that the separation 
of a few upward muon candidates was possible. 

We present the results of two independent analyses. One uses the approximation of the likelihood 
function by a F-function with an exponential tail [ 25 1 , the other the approximation by a Gamma 
function with an absorption term ||2^ (see section 6.3). 



Both analyses apply separation criteria which are obtained from a stepwise tightening of cuts on 
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different parameters, in a way which improves the signal-to-fake ratio given by the MC samples. 
Since the MC generated samples of downward-going muons (a few million events) run out of 
statistics after a reduction factor of about 10^, further tightening of cuts is performed without 
background-MC control until the experimental sample reaches the same magnitude as the MC 
predicted signal. 

For both analyses, the full experimental data set of 1996, starting with Feb. 19th and ending 
with Nov. 5th, was processed. It consists of 3.5 • 10^ events. 

Analysis 1 

In a first step, a fast pre- filter reduced this sample to a more manageable size. It consists of a 
number of cuts on quickly computable variables which either correlate with the muon angle, or 
which to a certain degree distinguish single muons from the downgoing multi-muon background 
events like, e.g. a cut on the zenith angle from a line fit |24|, cuts on time differences between OMs 



at different vertical positions, and topological cuts requesting a minimum vertical elongation of 
the event. 

These cuts reduce the size of the experimental data sample to 5.2%, the simulated atmospheric 
muons to 4.8% and simulated up-going events to 49.8%. 

Simulated up-going events and experimental data have been reduced by further cuts: 

• At least 2 strings have to be hit (this condition relaxes the standard condition ">3 strings" 
and increases the effective area in the vertical direction). 

• The events were reconstructed below horizon, i.e. 9 > 90". 

• log{Ctime)/Nhtt > -6. 

• a > 0.15 m/nsec, where a is obtained from a fit to 2:j = a • ij -|- /3 and Zi, ti being the 
z coordinates and times of the hit OMs. 



Fig. 2C shows the distribution of the number of direct hits, Naij.{15), of all events passing 



these cuts. The highest cut in Ndirect survived by any experimental event is Ndirect ^ 6. The 



two surviving events are shown in fig. 21. The Monte-Carlo expectation for upward muons from 
atmospheric neutrinos is 2.8 events, with an uncertainty of a factor 2, mostly due to uncertainties 
in the sensitivity of the detector after all cuts. 

Analysis 2 

The 3.5 • 10^ experimental events were compared to 3.5 • 10^ MC events from atmospheric down- 
going muons which correspond to 2 days effective live time. The MC data set for upward muons 
from atmospheric neutrino interactions ^^ consists of 2.5 • 10^ events triggering AMANDA-B4 
~ corresponding to 1.7 years effective live time. 

In order to separate neutrino induced upward muons, we applied a number of successively 
tightened cuts in the variables defined in section 6.4. This procedure reduced the experimental 
sample to the expected signal sample after the following cuts: 
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Figure 20: Number of events surviving pre-filter and additional cuts as a function of Ndir{i5). Solid 
line: 6-niontli experimental data, dashed line: 6-montli expectation from atmospheric neutrinos. 

1. reconstructed zenith angle 9 > 120°, 

2. speed of the line fit 0.15 < |i/| < 1 m/nsec, 

3. "time" likelihood log{Ctime) / {Nhit — 5) > — 10 (i.e. normalizing to the degrees of freedom 
instead of the the number of hit PMTs), 

4. "hit" likelihood logiCh^t)/iNhit - 5) > -8, 

5. number of direct hits for 25 usee window, Nf^ir{25) > 5, 

6. number of direct hits for 75 usee window, Ndir{7b) > 9, 

7. projected length of direct hits for 25 nsec window, LdirC^^) > 200 m, 

8. absolute value of the vertical coordinate of the center of gravity I^cogI < 90m 
(with the center of the detector defining the origin of the coordinate system) . 

Three events of the experimental sample passed these cuts, corresponding to a suppression factor 
of 8.9 • 10~^. The passing rate for MC upward moving muons from atmospheric neutrinos is 
1.3 % which corresponds to 4.0 events in 156 days. The corresponding enrichment factor is 
0.013/(8.9 • 10~^) « 1.5 • 10^. One of the three experimental events was identified also in the 
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Figure 21: The two experimental events reconstructed as upward muons. left: ID 8427997, right: ID 
4706870. The line with an arrow symbolizes the fitted muon track, the lines from this track to the OMs 
indicate light pathes. The amplitudes are proportional to the size of the OMs. The numbering of the 
OMs refers to the time order in which they are hit. 



search from the previous subsection. A second event with Ni^, 
analysis, with the exception of the N^i^ cut. 



5 passes all cuts of the previous 



In order to check how well the parameter distributions of the events agree with what one ex- 
pects for atmospheric neutrino interactions, and how well they are separated from the rest of 
the experimental data, we relaxed two cuts at a time (retaining the rest) and inspected the 
distribution in the two "free" variables. 
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Figure 22: Top - distribution in parameters Ldir{'25) vs. Ndir{75), bottom: distribution in the "com- 
bined" parameter S = Ndir{75) -Ldiri^^) / 20. The cuts applied to the event sample include all cuts 
with the exception of cuts 6 and 7. 



Fig. |2^ shows the distribution in Ldir{'^^) and Ndir{75). The three events passing all cuts are 
separated from the bulk of the data. At the bottom of fig. ^, the data are plotted versus a 
combined parameter, S = {Ndir{75)-2) ■Ldir{'25)/ 20. In this parameter, the data exhibit a nearly 
exponential decrease. Assuming the decrease of the background dominated events to continue 
at higher S values, one can calculate the probability that the separated events are fake events. 
The probability to observe one event at 5 > 70 is 15%, the probability to observe 3 events is 
only 6 • 10-"^. 
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Fig. 23 shows the distribution when \v\ and Lrfj.r(25) are relaxed. The 3 events are marked by 
arrows. There is one additional event at high LdirC^^), which, however, has a somewhat too 
small \v\. The 3 events fall into the region populated by MC generated atmospheric neutrino 



events passing the same cuts (bottom of fig. 23). We attribute the lack of experimental events 
between L[iij.{25) ~ 150-200 to statistical fiuctuations. 
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Figure 23: Distribution in parameters \v\ vs. Ldir{'25), after application of all cuts with the exception 
of cuts 2 and 7, which have been relaxed, top: experimental data, bottom: signal Monte Carlo sample 

Due to CPU limitation we could not check the agreement between experimental data and at- 
mospheric muon MC down to a 8.9 • 10~^ reduction. However, down to a reduction of 10^^, the 
disagreement does not exceed a factor of 3. A less conservative estimate of the accuracy of the 
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signal prediction can be obtained by replacing all dedicated cuts for 6 > 90° by the complemen- 
tary cuts for 9 < 90°. We observed a better-than-10% agreement between experimental data 
and MC after all cuts. In conclusion we estimate the uncertainty in the prediction of upward 
muon neutrinos to be about a factor 2. 

Table ^i^ summarizes the characteristics of the neutrino candidates identified in the two analyses. 
Table 1: Characteristics of the events reconstructed as up-going muons 



event ID -^ 


147 742 


4 706 879 


2 324 428 


8 427 905 


NoM 


13 


14 


15 


8 


^^ string 


3 


4 


3 


2 


\og{l./{Nhit - 5)) 


-8.3 


-8.5 


-8.0 


-11.2 


Orec, degrees 


168.7 


165.9 


166.7 


175.4 


(t)rec) degrees 


45.8 


274.2 


194.1 


- 



We conclude that tracks reconstructed as up-going are found at a rate consistent with that 
expected for atmospheric neutrinos. The three events found in the second analysis are well 
separated from background proving that, even with a detector as small as AMANDA-B4, neu- 
trino candidates can be separated within a limited zenith angle interval. Meanwhile, a few tens 
of clear neutrino events have been identified with the more powerful AMANDA-BIO telescope. 
They will be the subject of a forthcoming paper. 



10 Conclusions 

We have described the design, operation, calibration and selected results of the prototype neu- 
trino telescope AMANDA-B4 at the South Pole. 

The main results can be summarized as follows: 

• AMANDA-B4 consisting of 80 optical modules (-|- 6 OMs for technology tests) on 4 strings 
has been deployed at depths between 1.5 and 2.0 km in 1996. Seven of the OMs failed 
during refreezing. We have developed reliable drilling and instrumentation procedures 
allowing deployment of a 2 km deep string in less than a week. In the mean time the 
detector has been upgraded to 302 (AMANDA-B4, 1997) and 424 (1998) optical modules. 

• The ice properties between 1.5 and 2.0 km are superior to those at shallow depths. The 
absorption length is about 95 m and the effective scattering length about 24 m. 

• The original calibration accuracy reached for geometry and timing of AMANDA-B4 was 
about 2 m and 7 nsec, respectively. With the upgrade to 10 strings, these values have been 
improved to 0.5 - 1.0 m and 5 nsec. 

• We have developed proper methods for track reconstruction in a medium with non- 
negligible scattering. With tailored quality cuts, the remaining badly reconstructed tracks 
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can be removed. The quality of the reconstruction and the efficiency of the cuts improve 
considerably with increasing size of the array. 

• Geometry and tracking accuracy of AMANDA can be calibrated with surface air shower 
detectors. The mismatch between showers detected in the SPASE air shower array and 
muons detected with AMANDA is about 4 degrees, in agreement with Monte Carlo esti- 
mates of the angular accuracy. 

• The measured angular spectrum of the intensity of atmospheric muons is in good agreement 
with other experiments and with model calculations. 

• First neutrino candidates have been separated with AMANDA-B4. The identification of 
upward muon candidates with an array of only 73 operating 8-inch PMTs is a demonstra- 
tion that deep antarctic ice is an adequate medium for doing neutrino astronomy. 

Amanda-B4 is a first step towards a large neutrino telescope at the South Pole. A ten-string 
array, AMANDA-BIO, has been taking data since 1997. Presently, BIO data are analyzed, and 
tens of clear neutrino candidates have been extracted, with a threshold of typically 50 GeV. 
The construction of AMANDA-II, a 30 000 m^ array, is underway. The long-term goal of the 
collaboration is a cube kilometer detector, ICECUBE. 
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